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Abstract 

The relativistic corrections in the Hamiltonian for a particle in a uni- 
formly rotating frame are discussed. They are shown to be negligible 
in the case of ultra-cold neutrons (UCN) in the Earth's gravity. The 
effect, on the energy levels of UCN, of the main term due to the Earth's 
rotation, i.e. the angular-momentum term, is calculated. The energy 
shift is found proportional to the energy level itself. 

PACS codes: 

03.75.-b Matter waves 

03.65.Pm Relativistic wave equations 

04.25.Nx Post-Newtonian approximation; perturbation theory; related 
approximations 
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1 Introduction 



The Granit experiment has verified the quantization of the energy levels 
of ultra-cold neutrons (UCN) in the Earth's gravity field, in measuring the 
transmission of UCN through a horizontal slit [TJ. New, more precise exper- 
iments are being foreseen, in order to accurately determine the energy levels 
[2]. The aim of this paper is to calculate the main effects of the Earth's 
rotation on the gravitational energy levels of UCN. The effect of the Earth's 
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rotation on the quantum-mechanical phase of neutrons has been detected in 
the neutron interferometry experiments of Werner et al. [3] . Their theoret- 
ical interpretation is based on the non-relativistic Hamiltonian for a neutral 
particle in the inertial/gravitational field of the uniformly-rotating Earth: 

~2 

H = ^ (jj.L + V , LErAp, p = -iHV, (1) 

2m 

with Vo = —mU^, where Uj$(r) is the positive Newtonian potential, and 
where cj is the constant rotation velocity vector. The gravity acceleration is 
g = +V£/n, thus Vo ~ — mg.x in the locally-uniform Earth gravity, where 
x = r — i"o, with ro taken, say, on the lower side of the slit and at its entrance, 
in the Granit setting (which is described in Refs. [TJ HJ [3]). 

In Sect. |2j we shall take into account the relativistic effects which occur 
for a Dirac particle, in order to see if they are large enough to be detected. 
To this aim, we shall use mainly a study, by Hehl & Ni [6], of the post- 
Newtonian Hamiltonian of a Dirac particle in a purely inertial field, and also 
a previous work |7J on the energy levels of a Dirac particle in a purely grav- 
itational field. In the last part (Sect. [2D, we shall apply these results to the 
stationary states of UCN in the Granit experiment: we shall first estimate the 
relativistic effects, and then we shall calculate the main effect of the Earth's 
rotation, which is the (non-relativistic) effect of the angular-momentum term. 



2 Hamiltonian for a Dirac particle in an Earth- 
based laboratory 

2.1 Statement of the problem 

The quantum equation for an elementary relativistic particle with spin | 
in a flat spacetime is the Dirac equation. In the case with a gravitational 
field, which is described by a curved spacetime, the Dirac equation is usually 
modified to the form derived independently by Weyl and by Fock in 1929, 
hereafter the Dirac-Fock-Weyl (DFW) equation. Being a generally- covariant 
equation, the DFW equation may also be used to describe the effects of a 
non-inertial frame. However, for a neutral particle obeying Dirac's equation, 
the presence of an electromagnetic field has no effect whatsoever. This is 
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contrary to the observed existence of a magnetic moment for the neutron, 
detected in magnetic scattering experiments: these gave the first indication 
that the neutron might be non-elementary, as is of course admitted in the 
framework of the quark model. 

Nevertheless, we may admit as a working assumption that, in the ab- 
sence of an electromagnetic field, the inertial/gravitational effects on neu- 
trons are correctly described by the DFW equation. Our task then is to look 
for an approximate Hamiltonian, suitable to describe accurately enough the 
effects of a weak inertial/gravitational field on a slow particle. There has 
been some amount of work devoted to the subject (see e.g. [HI El QUI EH]), 
even though, as noted by Varju and Ryder [T2], "there is unfortunately some 
disagreement among [the different] papers." This may be partly explained 
by noting that the operator transformations of the Foldy-Wouthuysen type, 
which are used in Dirac theory to select the positive-energy states, are not 
uniquely defined [11] . In addition: i) the approximation scheme was not made 
explicit in these works, and ii) these works were concerned with the Hamil- 
tonian for the four- dimensional complex wave function, whereas it turns out 
that, as for the flat-spacetime Dirac equation, the solutions of the station- 
ary DFW equation have only two independent complex components [7\. In 
the asymptotic post-Newtonian scheme used in Ref. [TJ, the positive-energy 
states are automatically selected by assuming that any solution of the 
eigenvalue problem for the stationary energy levels should have an expansion 

EW=mc 2 + 0(A), (2) 

where A -C 1 is the gravitational field-strength. In that way, one avoids the 
ambiguity which may occur in the Foldy-Wouthuysen-type transformations. 
This study leads to an explicit equation for the stationary energy levels in 
a gravitational field, which shows that, for UCN, the corrections brought by 
the DFW equation to the Schrodinger equation in the potential mgz are ut- 
terly negligible |7J. This result has been since confirmed by an independent 
study [13]. 

However, these two studies [13] were based on the assumption of a 
static metric (as were the works [HI QUI EEE])- Hence, they do not involve 
the effect of the Earth's rotation, because the transition to a rotating frame 
implies a non-diagonal metric, goj ^ 0. One might try to adapt the foregoing 
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work [7J to a such "stationary metric." Our aim here is to estimate the "rel- 
ativistic" corrections, i.e., the difference between the predictions got either 
with the DFW equation or with the non-relativistic Schrddinger equation 
based on the Hamiltonian ([1]). Anticipating the smallness of these correc- 
tions, we may consider separately i) the "purely gravitational" corrections, 
i.e., those corrections which would be there if the Earth's gravitational field 
was static — the foregoing studies [T3] showed that these corrections are 
negligible for UCN — and ii) the "purely inertial corrections", which would 
be there if the rotating Earth had no gravitational field. The latter ones 
can be obtained from the work by Hehl and Ni [6], which is devoted to the 
inertial effects in a flat spacetime. 

2.2 Post-Newtonian Hamiltonian in a non- inertial frame 
in flat spacetime 

The work [6] is based on considering a (point-like) observer subjected to 
non-inertial (accelerated) motion in flat spacetime, thus in the absence of 
a gravitational field; due to the assumed rotation of an orthonormal triad 
(e^) that the observer carries with him, this approach also takes into account 
the rotation of the local reference frame. According to Kiefer and Weber 
[H], this approach is able to describe also the local effect of the gravitational 
field, by substituting — g for the three-acceleration a. Hehl & Ni [6] de- 
fine several coordinate systems and tetrads, one of which, (x^) [following the 
rotation of the triad (ej) and the corresponding tetrad (e M )] being qualified 
thus: "Such a local coordinate system is what we actually use in our labora- 
tory." Writing the DFW equation in that system, and using three successive 
Foldy-Wouthuysen transformations, they get the following form of the DFW 

1 These authors write: "by replacing a with g." Of course, the correct sign depends on 
which acceleration one refers to. Recall that, in a non-inertial frame F in non-relativistic 
mechanics, the inertial force on a particle is Fi = m(a — ao), where a is the acceleration of 
the particle in F and ao is its acceleration in an inertial frame — allowing to write Newton's 
second law in F as F = Fi + Fo = ma, with Fo the force seen in an inertial frame. A 
particle that remains at rest in F (relative velocity u = 0, and a = 0) is hence subjected to 
the inertial force Fj = — mao(r, u = 0). Thus, the absolute acceleration field ao(r, u = 0) 
of the observers bound to F plays the role of the opposite of a gravitational attraction 
field, g = -ao. 
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equation [6] (writing (3 = 7 , the standard-representation Dirac matrix): 

ih^ = HpwV' , (3) 



with 

2m 1 2m c 2 



H FW = /3mc 2 - 7— A + /?m(a.r) - w.L - w.S - 7— V— .V 



zft 2 

-cr. (a A V) + higher — order terms, (4) 

4mc 2 

with S = |ftcr, and where cr denotes the space "vector" made with the Pauli 
matrices cr J ' — or rather, with the matrices 

*-(#)• (6) 

since here the wave function ip has four complex components. In Ref. [T4], an 
additional term — 8 ^ c2 A 2 is present, and the remaining higher-order terms 
are stated to be 0(c~ 3 ). Moreover, in the last term above, cr.(a A V) is 
replaced by Per. (a. A V), which has the same effect on the "positive-energy 
spinor" (p, in a decomposition ip = (0, x) of the four-components spinor ip. 
Indeed, we have 

/3V =(&-*)■ (6) 

Thus, the second, third and fourth terms in the post-Newtonian DFW Hamil- 
tonian correspond exactly to the non-relativistic Hamiltonian ((Tj), if we 
put a = — g. This must be the case, since the other terms in (jll) are either 
the spin effect — u>.S, necessarily absent with the scalar Hamiltonian (CQ), 
or the relativistic effects. However, it is not necessary for us to mimic the 
gravitational field in a flat spacetime, since in this subsection we want to 
describe specifically the relativistic corrections to the non-inertial effects in a 
flat spacetime. Thus, what we actually should recover in the non-relativistic 
limit is the non-relativistic Hamiltonian ([TJ with V = 0, corresponding to 
a = in (jl]). It means that, for the application of the work [6], the "non- 
inertial observer" is on the axis of rotation of the Earth, and his "local triad" 
is just the one that rotates as does the Earth. 
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3 Application to the stationary states of UCN 
in the Earth's gravity 

3.1 Estimate of the relativistic effects 

Searching for stationary solutions ip(t, r) = <fi(t) A(r) of the Schrodinger-like 
equation (j3J), one gets energy states: 

R FW A — E A. (7) 

Setting 

A(r) = (^(r), X (r)), (8) 
we have explicitly from (HI) , for the independent 2-spinor ip: 

h 2 h 2 ( i 
H FW w = mc 2 A + m(a.r) — w.(L + S) a.rA + a.V + -<r.(a A V) 1 p. 

2m 2mc 2 \ 2 

(9) 

Let us briefly compare Eq. based on the work [6J, to other recent ap- 
proaches to the positive-energy PN Hamiltonian for a Dirac particle. Thus, 
we may compare ([9} with Eq. (52) of Boulanger et al. [13] , also based on a 
Foldy-Wouthuysen transformation, by assuming u> = here (since the rota- 
tion is not accounted for in Ref. [13]), and by cancelling the magnetic field 
B and assuming a = 1, b = in the latter work — so that Eq. (52) of Ref. 
[T3] corresponds to the Rindler metric 

ds 2 =(l + i^y c 2 dt 2 - dx 2 - dy 2 - dC,\ (10) 

relevant to a uniformly accelerated system. We find then an exact agreement 
if we substitute g = (0, 0, g) for a here, accounting for the fact that S = her/2, 
with h = 1 in Ref. [13]. This was expected, for the metric (jTOl) coincides 
at PN order with the metric used by Hehl & Ni [6], when uj = 0. To get a 
description of the rotation effects using the approach of Ref. [T3], one could 
substitute the Kerr metric for the Rindler metric. We may also compare 
with the positive-energy PN Hamiltonian that appears implicitly in Eq. (67) 
of Ref. [7], which corresponds to the metric 

ds 2 =(l- 2 c 2 dt 2 -fl + 2 (dx 2 + dy 2 + dz 2 ). (11) 



c 2 J \ c 2 
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In Ref. [7J, the positive-energy states are automatically selected by the 
Ansatz ([2]), as already noted. Naturally, the non-relativistic part in (Q coin- 
cides with that in Ref. [7J, if we put uj = and substitute = g.r for a.r, 
with g = —a. (Cf. Note 1.) With this substitution, the two last (relativistic) 
terms in are identical to corresponding terms in Eq. (67) of Ref. \J\ , up 
to positive factors, but the term with a.rA is missing, and a relativistic term 
involving A[/"n has to be added. Of course, in that case, an exact agreement 
was not expected. 

Let us now apply this to UCN on the rotating Earth. In Eq. a is 
then the absolute acceleration (denoted a in Note 1) of the non-inertial ob- 
server, which is at rest on the axis of the rotating frame. Thus a = in fact, 
as noted at the end of the foregoing Section. That is, if one computes the 
energy levels E m = E — mc 2 for the merely-rotating Earth, i.e., neglecting its 
gravitational field, the mere correction brought by using the PN Hamiltonian 
([HD with a = instead of the non-relativistic Hamiltonian ([T]) with Vq = 0, is 
the spin-rotation coupling term — oj.S. This is not really a relativistic effect, 
but one that arises due to the fact that the scalar Hamiltonian ([1]) does not 
involve the intrinsic spin. In other words, the relativistic effects of the pure 
rotation are simply zero at the PN approximation. On the other hand, recall 
that the relativistic effects of the pure gravitational field of the Earth (i.e., 
neglecting the rotation, but accounting for the global gravitational field) on 
the energy levels of UCN have been independently assessed [7J [13] and found 
fully negligible: they amount to less than ~ 3 x 10 -53 J. Our net result is 
thus that all relativistic effects are negligible for UCN in the Earth's gravity. 

The effects of the non-relativistic gravitational potential — mg.x are 
described in detail in Refs. [H HJ E], which, however, neglect the Earth's ro- 
tation. Hence, apart from the negligible relativistic effects, the perturbation 
to the stationary states as described in Refs. [TJ HI [5] is due only to the terms 
—cj.Li and — oj.S in the Hamiltonian which enters the eigenstate equation 
((7j). We have uo ~ 7.27 x 10 _5 s _1 . The spin-rotation coupling term — uj.S 
(first noted by Mashhoon [15]) is hence approximately 

| uj.S 1(T 39 J, (12) 

which is still small, though less small, as compared with the non-relativistic 
energies E nT , that start from 1.4 peV |TJ, which is ca. 2 x 10 -31 J. However, 
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it may be argued |1 6J that this coupling term has been indirectly measured 
in an experiment involving nuclear spin precession of atomic mercury |17j . 



3.2 Effect of the angular-momentum term on the en- 
ergy levels 

Let us now estimate the effect of the angular-momentum term <5H = — uj.L, 
which is the non-relativistic effect of the rotation, and is called "Sagnac 
effect" [3J EH [15] • Taking cylindrical coordinates p,9,( with the axis 0( 
being the rotation axis, thus parallel to vector uj, we have 

— uj.L(p = ihw.(r A VV) (f3) 

and 

u = ue r = pe p + (e V<p = —e p + -— e e + — e f , (14) 

whence 

(SB)(p = —u>.Lip = ihijj-— = ihu}p(V<p).e e . (15) 

do 

To assess this, we take for (p the unperturbed non-relativistic stationary state 
in the Earth's gravity [T], and we adopt Cartesian coordinates x,y,z from 
point A, with radius vector r , taken on the lower side of the slit [see after 
Eq. P])], the axis Ax being the West-East direction: e x = e 9 , and the axis 
Az being the local vertical. The stationary state is assumed to have the form 

m 

y>(x) = <p v (z)e i{klX+k2y \ (16) 

hence from (1151): 



(5H)to = iftwp— = —hukipip = —umuipip, (17) 
ox 

with 

Mi = u.e$, u = h{kie x + A^e^/m. (18) 

We note that the effect of the perturbing Hamiltonian <5H = — u>.L cannot 
be reduced to that of a mere additional potential SV, since, even for the 
particular wave functions f|T6l) . the ratio ((5l£)(p)/(p depends on cp through 
U\ = hki/m. We note also that this relevant velocity U\ is the West-East 
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velocity. (Compare Ref. [15].) 



In the neighborhood of the reference point A, the dependence of p (the 
distance to the rotation axis) on the East- West coordinate x is second-order. 
Calculating thus p in the "North-vertical" plane Ayz, we define the angle (3 
between Az and the radius vector AM. This adds to the latitude angle a 
to make the angle between the parallel to the equator and AM in the plane 
Ayz. We have cos/? = z/(y 2 + z 2 ) 1/2 , sin/? = y/(y 2 + z 2 ) l/2 , hence we get (R 
is the Earth's radius and po = -Rcosa;): 

P — Po + (y 2 + z 2 ) 1 / 2 cos(a + 0) = po + z cos a — y sin a. (19) 

We compute the modification (fT7|) of the energy level as a first-order pertur- 
bation: 

(5E) TOt ~ (<p | (5H)ip) = —uoraui J (p + z cos a — y sin a) \(p\ 2 dV, (20) 
with 



(po + 2 cos a — y sin a) |<p| 2 dV = po — sin a j y |<p| 2 dV" + cos a y 2;|<p| 2 dK 

(21) 

Due to the form of the wave function, the first integral (with y) does not 
depend on the quantum state. Thus, in (1211) . the contribution of the two 
first terms to 5E (!20l) is independent of the quantum state, hence may be 
omitted, because the energy levels are defined only up to a constant. That 
is, only the last term in (|2~T1) contributes: 

(8E) Iot ~ — umui cos a J z\(p\ 2 dV. (22) 



The integral in (1221) is the average value < z > n of z for the stationary 
state ip = <f n . If we use reduced quantities, 



z 



E { _ ( H 2 y _ fh 2 mg 2 



« S T- As 7' ' s feJ ' e n^J ' (23) 

we have exactly {[IE], eqs. (4b) and (18b)} 

< t >n= \\n- (24) 
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By (J23D, it follows that 



2 2 E 2 

< z >„= ~ZA n = --^ = -H n , (25) 

3 3 mg 3 

where H n = E n /{mg) is the classical turning point [5]. By (|22|) . this implies 
that the relative modification of the energy level E n is independent of E n : 

VKU ^ _ cos a 1Q _ 5 

(for Mi ~ +10 m/s and since uj ~ 7.27 x 10~ 5 rad/s, cos a ~ 0.71 and g ~ 10 
m/s 2 ). Hence, the effect of the Earth's rotation is not constant (as one finds 
by neglecting p — p ~ 10 cm with respect to p = -Rcosa ~ 5 x 10 8 cm, 
which seems justified at first sight), instead it is more important for high 
energy levels. The classical energy of the particle in the rotating frame of 
the Earth is [20]: 

-^classical = jmu 2 + V Q - ^m(uj A x) 2 = ^mu 2 + mgz - ^mu 2 p 2 . (27) 

Thus, this result may be intuitively understood in saying that, when the 
altitude z of the particle is increased, its distance p to the rotation axis is 
increased proportionally: 5p = 5z cos a < /) from Eq. (fl9l) . 



Let us estimate the state number n and the corresponding energy E n 
and height H n , for which the modification (5E n ) rot of the energy level due 
to the Earth's rotation becomes as large as the difference AE n = E n — E n _i 
between two successive energy levels. The WKB approximation, which is 
better at large n, gives E n oc x 2 / 3 with x n = n — | [HE], hence 

AE n oc (l - + 0(xnh ^ Ax„* , (5E n ) Iot ~ Bxl, (28) 



2 This was guessed by V. V. Nesvizhevsky (September 2006). His guess led me to do the 
calculation with a variable p, Eqs. ([19 ]) — (f2T |) . giving the result (|26|) . Moreover, one referee 
of the present paper suggested that "it might be helpful to provide an interpretation of 
the result," and the other one wrote that "it is not a surprise that the modification of 
the energy level due to rotation increase with the energy, it just reflects the fact that 
centrifugal force increases with distance" — thus suggesting the explanation that follows. 
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thus AE n ~ \(5E n ) mt \ when x n ~ A/ \B\. From E 7 ~ 6.044 and £ 6 ~ 5.431 
peV [5], we find A ~ 1.158 peV and, using also (125]) . 5 ~ —5.75 x 10~ 5 peV. 
Hence 

A£ n ~ |(5£ n ) rot | when n ~ 20139 ~ 20000, (29) 
corresponding to (c/. Table 1 in Ref. [5]) 

H n ~H 7 (—J ~ (58.945 x 207.25) /an ~ 12 mm, (30) 



or 



\ 2/3 



£ n ~ £ 7 f ^ ~ (6.044 x 207.25) peV ~ 1250 peV. (31) 

In the neutron trap, neutrons with different values and orientations of the 
horizontal velocity, thus with different values of the West-East velocity U\ 
entering Eq. (1221) . shall be present. Hence, for heights h of the horizontal slit 
of the order of 12 mm and larger, the higher-energy part of the corresponding 
spectrum, with H n ~ h, would be transformed to a continuum due to the 
Earth's rotation. 
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